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We analyze the effect of exchanging quasiholes described by Gaffnian quantum Hall trial state 
wave functions. This exchange is carried out via adiabatic transport using the recently developed 
coherent state Ansatz. We argue that our Ansatz is justified if the Gaffnian parent Hamiltonian has 
a charge gap, even though it is gapless to neutral excitations, and may therefore properly describe 
the adiabatic transport of Gaffnian quasiholes. For nonunitary states such as the Gaffnian, the result 
of adiabatic transport cannot agree with the monodromies of the conformal block wave functions, 
and may or may not lead to well-defined anyon statistics. Using the coherent state Ansatz, we find 
two unitary solutions for the statistics, one of which agrees with the statistics of the non-Abelian 
spin-singlet state by Ardonne and Schoutens. 



I. INTRODUCTION 



The study of fractional quantum Hall (FQH) liquids 
has lead to the discovery of rich classes of topological 
phases. The rapid progress in the theoretical exploration 
of this particular niche of correlated electron systems 
was made possible by the development of a set of princi- 
ples for the construction of "fixed point" or "prototype" 
wave functions for possible FQH phases, as pioneered 
by LaughlirP and then greatly expanded by other semi- 
nal contributions!^' This progress was compounded by 
the discovery of a link between these special wave func- 
tions and conformal field theory (CFT) P Via this link 
the topological quantum field theory describing a state 
can be inferred, which for some cases of interest is non- 
AbelianP The possibility of anyonic quasiparticles with 
non-Abelian braiding statistics may be the most spectac- 
ular implication of this field theoretic mapping, and its 
experimental verification remains a profound challenge to 
date, with topological quantum computing being a pro- 
posed application!^ Even in theory, a direct verification 
of the implied field theoretic mapping is nontrivial, and 
involves the calculation of adiabatic transport of quasi- 
particles using full many-body wave functions. This is 
still difficult in general, but has been possible for Abelian 
states^ and recently for p + ip wave superfluids and 
Moore-Read-type quantum Hall statesPtHI The general 
interpretation of analytic trial wave functions as describ- 
ing topological phases is thus well accepted whenever 
these wave functions are related to the conformal blocks 
of a unitary rational CFT. In this case, the field theo- 
retic mapping described above implies that the statistics 
of quasihole-type excitations, as defined through adia- 
batic transport, is directly encode d in the monodromies 
of the associated conformal blocks P^l This means that 
for such conformal block wave functions, adiabatic trans- 
port is the same as analytic continuation in the quasihole 
coordinates, and Berry phase effects are trivial. This 
then allows one to obtain the quasihole statistics directly 
from wave functions, but without explicitly calculating 
the Berry connection, or Wilczek-Zee connection, that 



defines adiabatic transport! 13 * 14 ! Indeed, in may cases of 
interest, this remains a profound technical challenge. 

In addition to wave functions arising from unitary 
CFTs, there is also considerable interest in analytic trial 
states that are similarly related to nonunitary CFTsP^ 
The physical interpretation of such states remains much 
more subtle. Here, the field theoretic mapping employed 
in the unitary case does not lead to a topological quan- 
tum field theory that can serve as the low energy ef- 
fective theory of the state in question. In particular, 
the conformal block monodromies cannot be interpreted 
to describe adiabatic transport, as they do not result 
in unitary transformations on states. In contrast, adia- 
batic transport describes (a limit of) the time evolution 
governed by a Hermitian Hamiltonian, and is therefore 
always described by a unit ary tra nsformation. For these 
reasons, it has been arguecPESEEIltn] that states obtained 
from nonunitary CFTs describe gapless critical points 
within the phase diagram of quantum Hall states, es- 
pecially in those cases where a local parent Hamiltonian 
for the state exists. Examples of the latter kind include 
the Haldane-Rezayi (HR) stated and the state known as 
Gaffnian, which has surfaced in the literature as early 
1993 through a fixed many-body clustering property,^ 
and which has been thoroughly discussed and was pro- 
posed to be critical using a CFT construction.^ 

The question arises what hidden orders can be identi- 
fied in such nonunitary states, whether they be remnants 
of topological orders or orders of a different kind. Un- 
fortunately, there is currently no efficient and universally 
applicable method to test for the topologically ordered 19 
nature of a state directly through the study of ground 
state properties. Much progress along these lines has re- 
cently been made through the analysis of entanglement 
spectra which are directly related to edge spectra. It 
has been argued that the edge spectrum of the Gaffnian is 
inconsistent with that of any unitary CFT, and that this 
contradicts the existence of a gap in the bulk spectrum,^ 
which is required for a topological phase. In principle, 
topological orders can be detected through nonlocal or- 
der parameters'^ though it remains difficult to explicitly 
construct such objects for general non-Abelian topolog- 
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ical orders in microscopic quantum Hall wave functions. 
The situation is similar in quantum magnetism. There, 
nonlocal operators detecting a topological phase can be 
directly constructed for toy modeliP^ defined on highly 
constrained Hilbert spaces where a gauge structure is ex- 
plicit (see Ref. f24] for a general discussion). However, 
such order parameters generally remain elusive in mod- 
els where similar physics is emergent within the low en- 
ergy sector of a larger Hilbert space (e.g., Ref. 25). The 
situation is much simpler in one-dimensional systems ex- 
hibiting a Haldane or Luttinger liquid phase. The hidden 
orders of these phases can be probed through non local 
objects measuring squeezed particle configurations J2SE3 
and their origin is quite manifest, e.g., in certain lim- 
its of Luttinger liquids where the wave function assumes 
a special factorized forml^SHSH For topological orders, on 
the other hand, the most direct probe that can, in princi- 
ple, be implemented at a microscopic wave function level 
is the study of the braiding statistics of localized elemen- 
tary excitations. 

In this paper, we are interested in some formal proper- 
ties of the Gaffnian state. We ask the question whether 
the Gaffnian trial wave functions may define some uni- 
tary anyon model through the holonomy calculated along 
exchange paths. Somewhat more physically speaking, 
this corresponds the the adiabatic transport of trial state 
quasiholes in the presence of a finite size gap. In- 
deed, this question is mathematically well defined. The 
quasihole trial states can be characterized as the unique 
zero-energy eigen states (zero modes) of a local parent 
HamiltonianpSES] For given quasihole configuration, the 
associated conformal block wave functions define a finite- 
dimensional subspace, which can be interpreted as a fiber 
over a point in the quasihole configuration space. The 
question is thus whether the holomomy associated with 
exchange paths in this configuration space induces well- 
defined statistics. It is clear from the outset that if this is 
so, the holonomies must be quite different from the con- 
formal block monodromies, since these holonomies give 
rise to unitary transformations on fibers. Physically, this 
is clear from the fact that these holonomies describe the 
adiabatic transport of quasiholes protected by a finite 
size gap. Mathematically, it follows from the fact that 
the connection on our vector bundle is a Wilczek-Zee 
connection defined in terms of a physical scalar product. 

The question defined above can be rigorously ad- 
dressed only by calculating the Wilczek-Zee connection 
from the given analytic wave functions. Unfortunately, 
we do not know how to do this for the Gaffnian state. In- 
stead, we will use this question as a testbed for a recently 
developed coherent state methocP^H to calculate adia- 
batic transport of quasiho les ba sed on the "thin torus" 
(TT)PES or "dominance"Ii2Il!Hm patterns of the wave 
function. This method has been shown to be quite effi- 
cient for a number of interesting states based on unitary 
CFTs, but has thus far not been applied to the nonuni- 
tary case. Our motivation to clarify the applicability of 
this method to a nonunitary state is twofold. A neg- 



ative result (no consistent anyon model) would further 
strengthen the case that the TT limit contains informa- 
tion about the gapped or gapless nature of the underlying 
state. This has been explored by one of us recently for the 
HR statepS though not with regard to statistics. On the 
other hand, if a consistent anyon model is obtained, we 
can argue that this is at least a very plausible scenario for 
the holonomies defined by the Gaffnian quasihole states, 
as we will further elaborate below. 

The coherent state method is based on the assumption 
of adiabatic continuity between states defined on a torus 
with arbitrary aspect ratio and corresponding states in 
the thin torus limit. It further rests on general assump- 
tions about a coherent state Ansatz for localized quasi 
holes in terms of adiabatically continued TT states. De- 
tailed arguments in favor of this Ansatz have been given 
in Ref. [34] Some of these arguments, in particular the 
justification for the factorized form of the Ansatz (see 
Eq. ([3]) below), also rest on a notion of locality, which is 
more subtle in a gapless state. We argue however, that 
the necessary assumptions still apply, as long as there is 
a finite gap in the charge sector of the system, indepen- 
dent of the existence of gapless neutral excitations. The 
scaling of the charge gap of the Gaffnian state has been 
discussed in some detail in Ref. 0H1 but at the moment, 
the question whether it remains finite in the thermody- 
namic limit has not been conclusively resolved to the best 
of our knowledge. 

There is a close connection between TT patterns and 
CFT fusion rules, which has been elaborated in Refs. |3HI 
and [501 In terms of the data used to construct an anyon 
model, there is thus some similarity between the present 
method, and the procedure of constructing _F-matrices 
consistent with given fusi on ru les, using the axioms of 
modular tensor categories! 51 * 52 ! Indeed, the two methods 
have so far given consistent, if not always identical, re- 
sults in the unitary casesP^ Some differences between 
these two approaches are worth noting. While the con- 
nection with tensor category theory can be physically jus- 
tified using the general framework of local quantum field 
theory,^ the assumptions used in the coherent state ap- 
proach are not field theoretic in character. Indeed, the 
i^-matrix, although it could be ultimately constructed, 
does not directly appear in this approach, and none of 
the defining consistency equations of this approach make 
reference to it. Moreover, the coherent state method, be- 
ing ultimately based on adiabatic transport, could not 
in principle lead to a nonunitary anyon model in its cur- 
rent formulation. A secondary, but non-negligible moti- 
vation for this work is thus to shine further light onto 
the connection between these two different methods in 
the context of a nonunitary state. 

The paper is organized as follows. In Sec. |H] we will 
summarize our calculation, which leads to two closely re- 
lated anyon models of Fibonacci type. All details are 
relegated to two Appendices. In Sec. |III[ we discuss 
our results, and make contact with the non-Abelian spin 
singlet (NASS) of Ref. 54. In a future version of this 
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preprint soon to follow, we will also present new numer- 
ical studies on the zero mode spectrum of the Gaffnian 
state and its connection to root partitions. 



II. BRAIDING STATISTICS THROUGH THE 
COHERENT STATE METHOD 

Our method to extract the statistics of a state from 
thin torus patterns, based on the assumptions discussed 
in the Sec. |TJ has been documented in detail in a recent 
paperPS Here we will describe its setup for the Gaffnian 
state, and focus on differences that arise compared to the 
discussion of the [k = 3) Read-Rezayi state^ as given in 
Ref. M 

The thin torus patterns of the bosonic v = 2/3 
Gaffnian, and their relation to the underlying mini- 
mal model CFT, have been thoroughly discussed by 
ArdonneP^ These patterns can be identified in the usual 
wa jj32j|44] w ] aen wave functions on the torus are consid- 
ered in the formal limit of an "extreme aspect ratio", 
e.g. L x ^S> l,L y <C 1. (Here and in the following, we set 
the magnetic length £b equal to 1.) The patterns emerg- 
ing in this limit are occupancy numbers in a suitably 
chosen Landau level basis on the torus (see, e.g.,[M]), de- 
scribing the trivial non-interacting product state result- 
ing in this limit. For the six degenerate Gaffnian ground 
states, these patterns read 200200 . . . , 020020 . . . , and 
002002 . . . , which we will call "(200)-type" , or 011011 ... , 
101101 . . . , and 110110 . . . , which we refer to as "(011)- 
type". 

As usualpHSI we w ju assume that states defined on 
a "bulk torus" — L x ^> l,L y >• 1 — can be evolved 
adiabatically into the thin torus limit described above, 
where the evolution is governed by the Gaffnian par- 
ent HamiltoniarP^ with slowly varying aspect ratio. 
Likewise, states with quasiholes on the bulk torus are 
assumed to evolve into elementary domain walls be- 
tween ground state patterns. There are three ele- 
mentary domain wall strings: 100 and 001, which oc- 
cur between (200)-type and (Oll)-type ground states, 
and 010, which occurs between two different (011)- 
type ground states. These domain walls may link var- 
ious different combinations of ground state patterns, 
thus forming charge 1/3 solitons. Representative ex- 
amples are given by the strings . . . 01101100200200 . . . 
and . . . 011011010110110 .... All charge 1/3 domain wall 
patterns occurring in the TT limit are locally given by 
translated and/or inverted versions of these two types of 
strings. This establishes a notion of "fusion rules" in the 
TT limit P 

The thin torus states are naturally labeled by the po- 
sitions of the domain walls together with a topological 
sector label — i.e., a label identifying the sequence of pat- 
terns between the domain walls. We use the notation 
|oi, . . . , a n ; c, a) for the TT states appearing in the limit 
L x ^> l,L y <C 1, and \a\, . . . , a n ; c, a) for the states 
appearing in the L x <C l,L y 3> 1 limit. The latter 
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TABLE I. c = thin torus patterns for a two-quasi-hole 
Gaffnian state, and the offset functions of the associated do- 
main walls. The elementary domain wall strings are in bold, 
and the orbital positions, 3rij, are underlined. Patterns for 
c = 1(— 1) can be obtained by shifting each occupancy number 
one orbital to the right (left), and the shift functions obtained 
using fj{c,a) = fj(a)+c. 



are simple product states in a "rotated" or dual Lan- 
dau level basia^jES and must be well distinguished from 
the states \a\, . . . , a n ; c, a). Here, the a, label the or- 
bital positions of the domain walls. For the topological 
sector labels, we follow the convention of [34] where a 
distinguishes classes of sectors that are not related by 
translation, and for given a, c = —1,0,1 distinguishes 
the three translationally related members of each class. 
Below we will also use translational properties to define 
a unique convention for how the c labels are to be as- 
signed. As usuaPS the domain wall positions can change 
only by multiples of a certain "stride" within each topo- 
logical sector (here, multiples of 3), and are thus of the 
general form cij = 3rij + fj(c, a). fj(c, a) is an offset fac- 
tor that depends on the sequence of domain walls (i.e., 
the topological sector and the domain wall in question). 
For symmetric domain walls, its value is uniquely deter- 
mined by symmetry. However, for asymmetric domain 
walls, a certain ambiguity exists a priori in how to define 
the domain wall position precisely with respect to the 
adjacent orbitals. This is accounted for by the shift — or 
asymmetry — parameter s. See 1341 and table [T] for details. 

By the assumption of adiabatic conti nuity, the "bare" 
domain wall states, \a\, . . . , a n \ c, a) and \ ax, . . . , a n ; c, a), 
give rise to two different complete sets of zero energy 
states (zero modes) of the Gaffnian parent Hamiltonian, 
for any aspect ratio of the torus. These two mutually 
dual zero-mode bases are obtained by adiabatically evolv- 
ing the respective bare states by means of a slow change 
in aspect ratio, and are denoted by \a%, . . . , a n ; c, a) and 
|<2i, . . . , a n ; c, a) respectively (with the dependence on as- 
pect ratio understood). These types of states generally 
describe quasiholes that are localized in x or y, respec- 
tively, at the domain wall positions, and are delocalized 
in the perpendicular direction (respectively y or x) . Their 
completeness within the zero-mode space implies that 
states with localized quasiholes (in both x and y) can 
be obtained by forming proper linear combinations. For 
these coherent states, with n quasiholes localized at com- 
plex positions hj — hj x + ihj y , the following Ansatz has 
been motivated in previous works (see Ref. [3H and refer- 
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ences therein): 

n 

, KOj) |oi, ...,a n ;c, a) 

a±<...<a n j=l 

(1) 

n 

\il>c,a({h})) = N' ^ Jj0 QJ -(/i J ,Ka J )|ai,...,a„;c,a) 

a!<...<a„ j'=l 

(2) 

with the Gaussian amplitude form factor 



a ,j(hj,a,j) = exp -(Ajj, + 8{a,j)/n)Kaj - j(h jx - KCijf 

(3) 

and its dual counterpart (j) a j(hj,aj) = 
<t>a,j(-ih>j,aj)\ K ^.fi. HerejK = 2n/L y , k = 2n/L x , 
and S(a,j) can be showrPS to be or 7r, taking on 
the same value for symmetry-related (translation or 
inversion) domain walls. Af and Af' are normalization 
factors. 

It must be emphasized that these coherent states are 
valid^ only for well-separated quasiholes: Eq. (Til) is valid 
when \hj X — h ix \ >• 1 for all and Eq. 15]) when 
\hj y — hi y \ ^> 1 for all i,j. When both conditions are 
satisfied, Eqs. ([!]) and ^ describe the same zero modes 
for given quasihole positions, and must be related by a 
linear transformation: 

= £ <.c'.^m)W7Wvi)- (i) 



L, 
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o 



o 
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FIG. 1. Two quasiholes, one stationary and one moving across 
the boundary, on a torus which is shown in a "repeated zone 
scheme". The stationary quasihole is at the bottom of the 
figure. Initially the moving quasihole (which starts at the 
dotted-line circle) is ordered second, and the system is in 
some configuration a. After that quasihole moves across the 
boundary it becomes the first, and the system is in configu- 
ration g x (o). 




FIG. 2. The configuration oi is the one in which the leftmost 
(first) quasihole is bottommost, the next to the right (second) 
is the next above, and so on. 



Here, a labels different "configurations" of the quasi- 
holes, i.e., components of the quasihole configuration 
space that can be connected without violating the con- 
ditions \hj x — h ix \ ^> 1, \hj y — hi y \ >■ 1. We define these 
configuration-labels with respect to quasihole configura- 
tions in the infinite plane. Toroidal periodicity leads to 
an equivalence between various triples (a, (c, a), (d , a')) 
that will be taken into account later. The transforma- 
tion described by v? cc , a i({h}) must also be unitary, 
since states in different topological sectors are generally 
orthogonal. 

A lot is known about the transition functions 
u cc' a fr° m the properties of the zero-mode ba- 



sis states — |ai, . . . , a n ; c, a) and \ai, ■ ■ ■ , a n \ c, a) — under 
magnetic translations alone Their dependence on h, 
(c, d), and (a, a') separates into the following factorized 
form: 



C,C' .Q,Q 



({/*}) = u({h})M c<c ,&, 



(5) 



with u({h}) and M c c > fully determined by translational 
symmetry: 



u({h}) = exp 



jxhjy 



(6) 



M = 



73 



e 27ri(i-l)/3 -2iriL/3 
g— 2niL/3 2 
„27ri/3 



e 2jri/3 
e 2jriL/3 
2niL/3 p2jri(£-l)/3 



(7) 



The ft,- independent product M CjC i£° a , may be identi- 
fied as the topological S'-matrix of the problem. We will 
now summarize how the missing information about a , 
can be obtained within this formalism, and subsequently 
this information can be used to gain knowledge about 
braiding. For details we refer the reader to Ref. [5H 

Further constraints on the matrix a , can be obtained 
by moving (via adiabatic transport) the positions of some 
of the n quasiholes across the boundaries of a rectan- 
gular coordinate chart on the torus, which we fix once 
and for all (see Fig. [T]). It is important to note that 
this can be done while maintaining \hj x — h ix \ 3> 1 and 
\hj y — hi y \ 3> 1, and this leads to new consistency con- 
ditions on the £^ / matrix. We begin in a configuration 
a — e.g., the one shown in Fig. [2] — and carry out this pro- 
cess in two steps. First, we consider the adiabatic trans- 
port of the upper right particle along a path that moves 
it further to the right. As soon as the particle crosses 
the right boundary of our coordinate chart, the particle 
reappears on the left boundary. Viewed in this coordinate 
chart, the quasihole configuration has changed from a to 
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a new configuration g x {o~). At the same time, with the 
particle that was formally rightmost now being leftmost, 
the topological sector label (the sequence of patterns ap- 
pearing in the L y — > limit) of the state has changed. If 
a was the label of the topological (super-)sector initially, 
the new label F(a) is easily found from the associated 
thin torus patterns™ Thus, while the transition function 
expressing the duality relation Q is defined in terms of 
£a a' right before crossing the boundary of the coordinate 

chart, it is defined in terms of C 9 ^} 17 } , right thereafter. 
Is is easy to work out a condition enforcing continuity 
across the boundary, relating these two matrices,^ 

x diag[e- 27r4/ ^( Q )/ 3 ], (8) 

where the argument of diag[. . . ] specifies the a-th diago- 
nal entry of a diagonal matrix. The matrix B is defined 



°a,F(a')- 



A similar process can be considered 



where the topmost quasihole of a configuration a moves 
up, crossing the upper boundary of the coordinate chart, 
becoming the bottommost quasihole in the resulting new 
configuration g y (a). The analogous continuity condition 
on the £ matrices reads, 



£ ff * (<7) = diag[e 



-27ri/j(a)/3- 



x dia,g{e~ 2niL / 3 ~ iLS( - a ' n ^ 3 ~ i / 3 '^ k s ^ a - k ~>]B. (9) 

The purpose of Eqs. (|8| and ^ is twofold: They connect 
£ matrices for different quasihole configurations, which 
will be essential in describing braiding processes, and 
they impose constraints on any given £ matrix. To see 
the latter, focus on the quasihole configuration shown in 
Fig. [2j which we will now refer to as 07 . It is easy to see 
that a 1 is invariant under the two moves described above 
performed in succession, i.e., g y (g x {o'i)) = 07. Hence, 
Eqs. Q and ^ together constrain the matrix elements 
of 

Yet another way in which the transition functions and 
hence the ^-matrices are related for different quasihole 
configurations is by mirror symmetry. In the presence of 
a constant magnetic background field, mirror symmetries 
exist only in conjunction with time-reversal symmetry, 
which we will leave understood. The operators t x and 
T y associated with mirror reflections across the x and y 
axis are therefore anti-linear operators. It is then sim- 
ple to find relations between the ^-matrices using these 
operations,^ 



(S T )- 1 diag[e 
x diag[( 



iL/3j2 k S(a,k) 



,2W3£ fc (l+<5(a,fc)/?0/fc(")l 



(10) 



(:9Ty(<r) — c Jj a g[ e 27Ti/3E fc (l+'5(«,'s)/^)/fe(a)]^)* 

xdiag[e iL/3 ^" s(a ' k) }B T . (11) 

In the above, g Tx (a) is the configuration that is the mirror 
image of a under t x , and similarly g Ty (o~). The matrix 
B T appearing above is defined by (B T ) a , a r — 0^,^(0'), 



h c 



O 
hi 



FIG. 3. Exchange path for two quasiholes at positions hi and 
hi+i. First, the quasihole at hi+i is dragged to h a - There the 
coherent state representation is changed from the original ba- 
sis to the dual basis using H CT . That quasihole is then dragged 
to hb, and the state is changed back to the original basis using 
H CT . The quasihole at hb, now the i-th, is moved to h c , then 
both quasiholes are moved to their final positions: The quasi- 
hole at hi goes to h i+ \ and the quasihole at h c goes to hi. 
Using the coherent state Ansatz, this gives rise to Eq. (12 1. 



where F T (a) is the (super-)sector resulting from an in- 
version of the patterns associated with the (super-)sector 
a. Operating t x and t v in succession on a system in the 
configuration 07 gives 07 again, i.e., g T (g Tx (c/)) = 07. 
In this way Eqs. (10 1 and (111 allow us to further con- 



strain the elements of (f 1 . 

To braid two quasiholes, we adiabatically transport 
them along the path in Fig. [3] To this end, we consider a 
configuration a where both quasiholes are adjacent both 
in terms of their x and their y coordinates, i.e., they 
can be exchanged without violating \h 



JX 



x\ > 1 or 

^> 1 for any other pair of quasiholes. The re- 
sult of the adiabatic transport can then be worked out 
directly from the coherent state Ansatz (111) given known 
transition functions, i.e., the matrices £°TI2}E1] it can be 
expressed as follows: 



( 1*0 \ 

1*2) 

V l*n) J 



i*AB2 CT (2 <T ')^ 



( 1*0 \ 

1*2) 

V l*n) J 



(12) 



Here, $ab is the Aharonov-Bohm phase, equal to the 
charge of a quasihole, —1/3, times the area enclosed by 
the braiding path. The quantities \^ a ) are the three- 
component column vectors with entires \ip c ,a)- ^ is the 
matrix £ a ® M (where we identify the states \i/j c ,a) with 
a formal tensor product basis |a) ® |c)). a is the initial 
configuration of the quasiholes, and a' the other configu- 
ration that occurs during braiding (Fig. [3]), i.e., the one 
obtained from a by crossing the line hj x = hi X or the line 

It turns out that the result of braiding is always block 
diagonal in the c labels, i.e., the braid matrix E! <T (S IT )t 
is of the form Xi( n ) ® Ic max xc max , where we call Xi( n ) = 
? <T (£ cr ) the "reduced" braid matrix associated with a 
counter-clockwise exchange of the i-th and i + 1-st of n 
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quasiholes. This fact is a direct consequence of trans- 
lational invariance. Moreover, Xi( n ) is found to be in- 
dependent of the initial configuration cr, as one would 
expect. 

The various constraint equations discussed so far still 
admit many solutions for the braid matrix. This is chiefly 
due to the fact that the asymmetry parameter s intro- 
duced above is still undefined. A final set of constraint 
equations comes from the imposition of certain locality 
constraints on the braid matrix that can b e motivated di- 
rectly from the thin torus patterns ) 33 l 34 l Said succinctly, 
we mean by locality that the result of braiding should 
only depend on the sequence of three ground-state pat- 
terns forming the two domain walls associated with the 
braided quasiholes, and that only the pattern in the mid- 
dle may change as a result of braiding. The constraint 
equations following from this, combined with the above 
symmetries, then lead to a discrete set of (usually inti- 
mately related) solutions for the statistics. We will dis- 
cuss the full set of constraint equations and their solu- 
tions in App. [X]for the two-quasihole case, and in App. 
[B] for the three-quasihole case. Here we will summarize 
the results of this calculation by giving the (reduced) 
braid matrices obtained from it for both two and three 
particles: 

Xl(2)= ^)(^))t 



e -2i7ra ( ^-l/2 



e -2i7ra^-l/2 

v- 1 



(13) 



Xi(3)=^ W3 (C (2 ' 1 ' 3) ) t 



e -2i7ra ( ^-l/2 



(14) 



where a = ±1/5, and tp is the golden ratio, (p — 
(1 + y5)/2. Here, the rows and columns refer to the 
a-supersectors of Tables [I] and [H] The parameter s is 
found to have one of two values: s = 2 — 3a/2. Just as 
in the case of the k — 3 Read-Rezayi state discussed in 
Ref. [Ml there are two solutions, one for each sign of a, 
which are related by an Abelian phase and complex con- 
jugation. In the above, we have also fixed a gauge degree 
of freedom associated with unitary transformations. 

Together with the locality constraint described above, 
these two matrices imply the result of braiding any ad- 
jacent pair in a state of n quasiholes. A "tensor repre- 
sentation" of the statistics just as discussed in Ref. [55] 
can then immediately be constructed in complete analogy 
with Ref. mi 



a 


Thin torus pattern 




Ha) 




1 


002OOU011011O1O11011011OO2002 


— s 





s 


2 


110110101101100200200200110110 


1 


-2 + s 


1 - s 


3 


110110020020020011011010110110 


-1 + s 


2 - s 


-1 


4 


110110101101101011011010110110 


1 





-1 



TABLE II. c = thin torus patterns for a three-quasihole 
Gaffnian state, and the offset functions of the associated do- 
main walls. The elementary domain wall strings are in bold, 
and the orbital positions, 3nj, are underlined. Patterns for 
c = 1(— 1) can be obtained by shifting each occupancy number 
one orbital to the right (left), and the shift functions obtained 
using fj(c,a) = fj(a)+c. 



III. DISCUSSION 

The two solutions obtained in the preceding Sec. [n] 
are related to one another simply by complex conjuga- 
tion and an overall Abelian phase. They describe Fi- 
bonacci anyons and are thus closely related to those ob- 
tained for the k = 3 Read-Rezayi (RR) state using the 
same mcthodP In essence, the solutions obtained from 
the RR patterns and those obtained here are the same 
up to an Abelian phase. This statement excludes "global 
exchange paths" on the torus involving processes such 
as the one depicted in Fig. [TJ as we will further explain 
below. This close correspondence between the solutions 
found from RR and Gaffnian patterns is a manifestation 
of level-rank duality between the associated SU(2) 3 and 
SU(3)2 fusion rules, respectively. This duality was also 
discussed by Ardonne^l in terms of domain walls. It is 
manifest in the Bratteli diagrams of Fig. |4j which in the 
present context represent the rules for domain wall for- 
mation between ground state patterns for the respective 
states. Note, however, that there is a three to two corre- 
spondence between the topological sectors in both cases, 
rather than one-to-one. This is so since each sector, i.e., 
each path in the Bratteli diagram, is threefold degenerate 
under translations in the Gaffnian case, but only twofold 
in the k — 3 RR case (at filling factor 3/2). The cor- 
respondence between Gaffnian and RR sectors is perfect 
if we limit ourselves to the n — 1 generators of the braid 
group cr^i+i, i= 1, . . . , n— 1, which exchange the i-th and 
i + 1-st quasihole. These generate the full braid group in 
the plane, but not on the torus. On the torus, these gen- 
erators leave certain subspaces of topological sectors in- 
variant, which all start and end in the same pattern in the 
topological sector label. These subspaces for the Gaffnian 
are then in correspondence with similar subspaces for the 
RR state, in the sense that there are isomorphisms be- 
tween them that commute with braiding, except for an 
overall Abelian phase. This correspondence, however, 
gets spoiled by the inclusion of the remaining generators 
on the torus, which mix the subspaces. This happens 
differently for the Gaffnian and the RR case, since there 
are six such subsectors in the former case, but only four 
in the latter. 
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n 



FIG. 4. Bratteli diagrams depicting the c = patterns in the 
RR state (top) and the Gaffnian (bottom). Valid topological 
sectors on the torus are represented by paths which start at 
the left, take one step right (either up or down) for each quasi- 
hole in the state, and end on the same type of ground state 
pattern — (200)-type or (Oll)-type — as they began, minding 
periodic boundary conditions. There is a one-to-one corre- 
spondence between the paths in the lower diagram and the 
patterns in Tables |I| and |H} an d also between the valid paths 
in the upper and lower diagrams. This latter correspondence 
is how the SU(2)3 and SU(3)2 rank-level duality manifests 
in terms of patterns. It should be noted that for each sector 
which corresponds to a path in one of these diagrams, there is 
for RR an additional sector related to the first by translation, 
and two additional translated sectors for the Gaffnian. 



The subtle differences between our solutions for the 
Gaffnian and the RR case on the torus are of a piece 
with the difference in overall Abelian phase. It is well 
known that the overall Abelian phase could in principle 
assume any value in planar geometry, but on the torus, 
it is constrained by the topological degeneracies charac- 
terizing the state. In the coherent state method, one 
source of phase differences is the factor i/3 in the co- 
herent state Ansatz, Eq. ([3]), which is generally related 
to the "stride" of the domain wall in a given topologic 
sector, which equals 3 in the present case and 2 in the 
RR case. This stride is of course identical to the center- 
of-mass degeneracy. In particular, one may see that the 
equations obtained from global processes such as the one 
shown in Fig. [T] are quite sensitive to this stride and the 
associated phase (see Eq. (All below). In view of the 
importance of these processes in our method, and the 
fact that they spoil the correspondence between Gaffnian 
and RR topological sectors as explained above, it may 
not be clear a priori that the consistency equations we 
obtain in both cases admit closely related solutions, in 
the sense discussed. That this is so can be traced back 
to the fact that the translational degree of freedom, c, 
decouples early on (see below Eq. (|5|), and the remain- 



ing a degree of freedom is fully analogous in both cases. 
This is how rank-level duality becomes manifest in the 
present formalism. For similar reasons, our solutions for 
the Gaffnian and RR states, which were both obtained 
at the maximum (bosonic) filling factor, could be gener- 
alized quite easily to lower filling factors. 

As emphasized initially, the Fibonacci-type solutions 
we obtained are distinct from the anyon model asso- 
ciated with the conformal block monodromies of the 
Gaffnian state. The latter describes so-called Yang- 
Lee anyons, whose relation to Fibonacci anyons and 
the associa ted G alois duality has enjoyed much inter- 
est recently? 56 * 57 ! Yang-Lee anyons are associated with 
nonunitary F-matrices consistent with the SU(3)2 fusion 
rules. It has been known for some time, however, that 
the same fusion rules admit unitary solutions of the Fi- 
bonacci type also; these are realized by the NASS state 
of Ref. [5U Indeed, it is not difficult to perform checks 
confirming that one of our solutions, that corresponding 
to s = 17/10, agrees exactly with the monodromies of 
the NASS state (including the overall phase). The dom- 
inance patterns of the NASS state have been discussed 
more recently,^ and it seems clear that the calculation 
presented here can be carried over to this state without 
essential changes. In all, this confirms once more for the 
case of SU(3)2 that the coherent state Ansatz produces a 
subset of all unitary anyon models consistent with given 
fusion rules. 

In the case of the unitary NASS state, our results 
support the usual conjecture that the holonomies associ- 
ated with adiabatic transport agree with conformal block 
monodromies. In the case of the Gaffnian, things are 
more subtle. We have argued that in this case, provided 
that the coherent state Ansatz is justified, the holonomies 
gives rise to a well defined Fibonacci anyon statistics, 
possibly (given the twofold ambiguity of our solution) 
identical to those of the NASS state. In particular, we 
believe that our Ansatz is indeed well justified if the gap 
of the Gaffnian parent Hamiltonian does not close in the 
charge sector. If so, this would presumably remain a for- 
mal property of the Gaffnian state that is not robust to 
perturbations, even ones that do not open a gap in the 
neutral sector. However, it might shine new light on the 
formal connections between the Gaffnian and the NASS 
state, which have been discussed previously.^ We leave 
more rigorous investigations into this matter for future 
studies. 



IV. CONCLUSION 

In this paper, we have analyzed the effect of exchang- 
ing quasiholes described by Gaffnian quantum Hall trial 
state wave functions via adiabatic transport, based on a 
coherent state Ansatz in terms of adiabatically continued 
thin-torus zero modes. We have argued that this Ansatz 
may correctly describe the holonomies of Gaffnian quasi- 
holes as long as the charge sector of the state remains 



gapped. In this case, we find that the statistics obtained 
are either closely related, or in fact identical, to those of 
the NASS state, though this property will not be topo- 
logically protected. This observation is of a kind with the 
statement that these states have the same fusion rules, 
though their conformal block monodromies are different, 
and are non-unitary in the Gaffnian case. Our results 
appear to confirm in yet another case that the coher- 
ent state Ansatz produces a set of unitary anyon models 
consistent with the fusion rules of the underlying state. 
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Appendix A: Two-quasihole Solution 

To work out the braid group for n quasiholes with 
general n, one needs to consider only braiding for pairs 
of quasiholes associated with all possible pairs of do- 
main walls, as given by all possible combinations of three 
ground state patterns. Locality then implies that all the 
other ground state patterns appearing in the topological 
sector label will not affect the result of braiding. To this 
end, we will solve for the reduced braid matrix in the 
simple cases of n = 2 (here) and n = 3 (in Appendix [b]) . 
Together, these results can be used to construct braid 
matrices for n-quasihole states, since these cases exhaust 
all possible sequences of three ground state patterns sep- 
arated by domain walls. 

The set of constraint equations to solve for the reduced 
braid matrix come from enforcing the unitarity of the 
£ matrices, the locality condition discussed in the main 
text, and the global path relations Eqs. ^ and We 
will begin by enforcing the latter. As discussed above, 
we can apply Eqs. Q and (|9| in succession to constrain 
£, ai because 07 = gy(gx{o'i)h Applying these two equa- 
tions, with the data from Table|lJ results in the constraint 
equation 



respectively. Equation (Al) is satisfied when 



CO 



( Ap" 1 \ / Ap 
Ap Ap- 1 

V -1/ V - 1 

(Al) 

where p = — exp [— 2m(l + s)/3], A = 
exp[-27ri(£/2+l)D/3], and D = or 1 if the S 
parameter for the 100-type domain walls is or n, 




(A2) 



Mirror symmetry, Eqs. (10) and (ITT]), can also produce a 
constraint equation because g T (g~ x (07)) = 07. However, 



in the case of two quasiholes, applying Eqs. ( 10 ) and ( 11 ) 
in succession results in the trivial equation, £ ai = . 

The process of solving for the reduced braid matrix is 
similar to the solution given for n = 2 in Ref. [Ml We 
gain the following equations by demanding that is 
unitary, 



= 2Ar,- D pCn^2 + V D P^ 



2 

3 > 



(A3a) 

e-/3 = 2?? V 3 i 2 - £ 33 2 , (A3b) 
= At,- V^ii 2 + Arr Vi 2 2 - A ?? VV, 

(A3c) 

= v~ D &i (VCn + AK12) + (A3d) 

where rj — exp (— 2ni/3). Two additional equations come 
from the requirement that braiding is local; as said above, 
this means that the result of braiding should only depend 
on the sequence of three ground-state patterns forming 
the two domain walls associated with the braided quasi- 
holes, and that only the pattern in the middle may change 
as a result of braiding. Imposing these locality consider- 
ations tells us that \i (2) must be of the form 



Xi(2) 



^■00 
• • 



(A4) 



where "-"s are unknown, potentially nonzero, matrix el- 
eme nts fo r which we will solve. By applying the form in 
Eq. (A4| to the matrix xi(2) = £ <Tl (£ <Tj ') t derived from 
adiabatic transport, the zero elements give two more in- 
dependent constraint equations, 



= rV(l+ ? V 1 i£i2-A77V 3 2 , 
= Vi (£11 - AK12) + fVi 3 £ 33 - 



(A5a) 
(A5b) 



Solving this system of six equations, (A3) and (A5), 



is formally similar to the solution in the Appendices of 
Ref. [3H so the details will not be repeated here. Just 
as in that reference, there are two solutions: A special 
solution in which p = ±i, 



f 2 - -n D e™/ 3 n- 1 



62 

• 2 
33 



£31 = 0, 



(A6a) 

(A6b) 
(A6c) 

(A6d) 
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which produces the reduced braid matrix 

/ 



Xl (2)=e 2 ™/ 3 




(A7) 



(but we will show in Appendix |B| that this solution is in- 
consistent with the equations from three-quasihole braid- 
ing), and the consistent solution, 



a ?A W3 
11 (1+P) 2 ' 
£12 = Afu, 

^^(p+jt 1 )^, 

which produces the reduced braid matrix 
Xi(2) = e—/ 3 x 



(A8a) 

(A8b) 
(A8c) 
(A8d) 
(A8e) 



/ p- 1 








P(P + P _1 -1) ±e' i7r£l / 3 (l -p)y / p + p- 

\ ± e - i7rI5 / 3 (i -p)VP+]F T p + p- 1 -! 

(A9) 



This two-quasihole reduced braid matrix is the same as 
that in Ref. [34jexcept for two features: The independent 
phase on the off-diagonal elements (though this will later 
be removed with a unitary transformation), and the over- 
all Abelian phase, which here is e _I7r / 3 and in Ref. IM1 
was e i7r / 2 . 



Appendix B: Three-quasihole Solution 

The solution for the reduced braid matrix of three 
quasiholes begins similarly to that for two quasiholes. 
We first constrain £ ai using global path relations, Eqs. 
^ and ([9]), and mirror symmetry, Eqs. ( |To| and (111. 
Applying the former two in succession and filling in the 
data from Table [TT| gives the constraint 



Ap- 1 




Ap 



\ ( 



V 



■1/ 



Ap" 1 



V 



Ap 



(Bl) 



can be applied in succession to £?' to give 



/I 



T] D 

ir D 



V 



\ 



1/ 



/I 



r)~ D 
if 



V 



(B2) 



1/ 



Equations ( Bl I and (B2| together constrain £ CT/ to be of 
the form 





61 


62 




64 








-AKl2 


-A-^14 


V 




A^12 




-AK14 


V 




-A-^41 


-Apfci 


£44 



The structure of this matrix is almost identical to the 
corresponding matrix £ 
D-dependent phases r\ L 



++ 



in Ref. [3U save that the 
and A are different. There, 



A was defined such that A 



A-rrD 



whereas here 



A 2 



= 2i7r_D/3 



i)~ ■ We might then find a braid ma- 
trix with a non-trivial dependence on D. However, the 
symmetry relations in Eqs. (J8|, (10 1, and (11) also 
have additional D-dcpendent phases compared to their 
iccrresponding forms in Ref. 1341 and we will see that 
by following the same steps as in that reference to find 
the braid matrix solution, all the D-dependent phases 
will conspire to cancel save for those on the off-diagonal 
elements which can be removed via a unitary transfor- 
mation. 

By requiring that £ <7 - r be unitary, we find the following 
constraint equations, 



l = |£ii| 2 + 2|£ 12 | 2 + |£i 4 | 2 , 



(B4a) 



= A& a e u * + A -1 p 2 £ii£ 12 * - p |£ 12 | 2 - p \£ u \ 2 , 

(B4b) 







i;i (tii -A : +i. u "in- 



1 = 3 ie 4 ii 2 + ie 44 | 2 • 



(B4c) 
(B4d) 



The locality of braiding tells us not only that some el- 
ements of Xi(3) must be zero, as was the case for two 
quasiholes, but also that the 2x2 block with off-diagonal 
elements must be equal to the equivalent 2x2 block in 
Xi(2). This is because the sequences of ground state pat- 
Vterns of the domain walls associated with the quasiholes 
;o be braided are the same for those two supersectors in 
;he two- and three-quasihole cases. In other words, Xi(3) 
nust be of the form 



-1/ 



where p is defined as in Appendix |aJ and A = 
exp [—2iri(L/2 + l)D/3], This definition for A is seem- 
ingly the same as that for A in Appendix [5] but in the 
case of two quasiholes L = 1 modulo 3, and here L = 
modulo 3. 



Xi(3) 



e -<f/3 x 



p(p + p 1 + 1) 

±e- inD / 3 (l-p) v /p + p~ 



p + p~ 



The two mirror symmetry equations, (10) and (11), 



P 
1 

(B5) 



P 
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if Eq. ( A9 ) is the correct reduced braid matrix for two- 



quasiholes, and of the form 



Xi(3) 



V 



p 

e-"/ 3 j 



(B6) 



with p — ±i, if Eq. (A7) is the correct matrix (which 
we will show is not the case). We find our constraint 
equations by equating the product Xi(3) = (£, ai )t, 
obtained from adiabatic transport, with the forms above. 
We can perform the two solutions in parallel by using 
only the elements of Eqs. (B5 1 and (B6 1 that are zero in 



both. This produces the constraint equations 



0=^x1 -A'Cp-lJpf^ 



2f 2 
V f 14 , 



= A(p-l)p 2 £ii£i2 + A 2 Ki 2 2 



= -U 



Vixx + A(p - l)£i2 



-T?P 2 

P ?14 > 
' £l4£44- 



(B7a) 
(B7b) 
(B7c) 



The set of constraint equations, (B4) and (B7|, is solved 
when 



J8i 



Cl1 " (1+P) 2 ' 
ei4 a = (P + P _1 -l)en a . 



which produces the reduced braid matrix 

Xl (3) = e -W3+i9i x 

/l 



a) 

(B8b) 
(B8c) 
(B8d) 
(B8e) 



additional D-dependent phase and the overall Abelian 
phase is different. 

We have yet to enforce consistency between the 2x2 
blocks in the two- and thrcc-quasihole braid matrices; 
to do so we equate Eq. (B9) to Eqs. (B5) and (B6) 
in turn. The latter produces a contradiction, because 
the off-diagonal elements of Eq. (B9) are not zero for 
p = ±i. Thus Eq. ( A6 1 is not a consistent solution for two 



quasiholes, and Eqs. (A7) and (B6) are not consistent 



braid matrices. Enforcing consistency between Eqs. (B5 ) 
and ( B9 ) implies that 



P + P 

where tp is the golden ratio, ip = (1 
words, 



p — exp ± 



(BIO) 
/5)/2. In other 



(Bll) 



or, if we define a = ±1/5, p = exp {iita). This is all the 
same as in Ref. 1341 However, here the s parameter is 
defined differently in terms of p than in Ref. 1341 we have 
defined p = — exp [— 2wi(l + s)/3], so the s parameter is 
also constrained to be 



2 _3a 
2 



(B12) 



Consistency also implies 



e 

2i6 



\ 



P 
: 1. 



(B13) 
(B14) 



p(l -p) 



Th e 



±e 



± e i02-irrD/3 p ^ p + p -l _ 1 



expressions for xi(2) and Xi(3) in Eqs. (13 1 and 



d/3 pVpT: 

e 2l02 (l—p) to ^qs. (A9) and ( |B9[ ) and have undergone a uni- 



respectively, have been simplified with respect 



(B9) 



Just as in Appendix [Aj Eq. ( B9 1 is the same reduced 
braid matrix as was found in Ref. [331 for three quasi- 
holes, except that here the off-diagonal elements have an 



tary transformation in which each state is multiplied by 
(-e i7r£, / 3 ) #200 , where #200 is the number of 200200 . . . 
strings in the thin torus pattern associated with that 
state. This unitary transformation removes the depen- 
dence on the unkonwn parameter D. 
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